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Abstract 

Inspired by experiments on the actin driven propulsion of micrometer sized beads we 
develop and study a minimal mechanical model of a two-dimensional network of stiff 
elastic filaments grown from the surface of a cylinder. Starting out from a discrete 
model of the network structure and of its microscopic mechanical behavior we derive a 
macroscopic constitutive law by homogenization techniques. We calculate the axisym- 
metric equilibrium state and study its linear stability depending on the microscopic 
mechanical properties. We find that thin networks are linearly stable, whereas thick 
networks are unstable. The critical thickness for the change in stability depends on the 
ratio of the microscopic elastic constants. The instability is induced by the increase 
in the compressive load on the inner network layers as the thickness of the network 
increases. 

The here employed homogenization approach combined with more elaborated mi- 
croscopic models can serve as a basis to study the evolution of polymerizing actin 
networks and the mechanism of actin driven motion. 

Keywords: A-microstructure, B-biological material, B-constitutive behavior, 
C-stability and bifurcation, homogenization 
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1. Introduction 

Dynamic filamentous polymer networks play an essential role in the mechanics of 
living cells. However, adequate constitutive equations, which take into account the 
growth history of the networks and the thereby developping prestresses, are difficult to 
obtain. One prominent example, which shall serve as a starting point for this paper, is 
the actin filament network. Actin filaments cross-linked by a special protein complex 
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Figure 1 : The actin network (shades in red) grows initially in a symmetric cloud (a) from a bead surface 
(black) until a spontaneous symmetry breaking occurs and the growth continues in an asymmetric fashion 
(b). Newer (older) network layers are in dark (light) red. 



called Arp2/3 comple)(Q play a key role in cell motility Jl]] as well as in the motion of 
cellular organelles J2] or pathogens QUI] inside their respective host cells. 

This actin driven motion has been reconstituted under minimal conditions [1] using 
beads Jsf], vesicles |@], and droplets 10]. The basic biochemical processes necessary 
for motion are by now well characterized (for a review see 0] and references therein). 
Briefly, under the action of several auxiliary proteins, actin polymerizes into a polar 
filament network, whereby filament elongation and filament linkage are restricted to a 
zone close to the surface of the propelled object, i.e. the internal interface (The filament 
network is highly porous and does not hinder the diffusion of monomers to the surface 
of the object.). Depolymerization takes place far away from the object, i.e. at the 
external interface. Therefore, new network material is inserted between the already 
existing network and the object and forces the older network away from the object 
[see Fig.Q](a)]. If this type of growth is occurring on curved surfaces, the network 
experiences large stresses and deformation as growth proceeds. In experiments with 
actin-propelled biomimetic objects (e.g. beads) O, la LZD it was shown, that initially 
the actin networks forms a symmetrical cloud around the object, until a spontaneous 
symmetry-breaking occurs [see Fig.[T|(b)]. Then the network is growing faster on one 
side than on the other and finally the object starts to move pushed by an actin comet. 
Thereby actin is continuously polymerizing in the contact region between the bead and 
the comet. 

Thus the (chemical) out-of-equilibrium process of polymerization/depolymerization 
(which is nevertheless stress dependent, see |@] and references therein) and the action 
of enzymes produce a network in a pre-stressed state, which does not relax on the 
timescale of symmetry-breaking and which depends on the history of growth 0, [Tol 
ll lL 1 1 2n . Several theoretical studies have addressed this problem either using discrete 
ll 5L 1 1 3ll or phenomenological continuous models lfl4j[l5 , 16 1 without, however, obtain- 
ing an adequate constitutive law for the network. A type of homogenization model has 



The Arp2/3 complex is a special protein complex which induces the nucleation of new actin filaments 
at already existing ones, thus triggering the formation of Y-junctions. It is activated by other enzymes (e.g. 
ActA, Wasp) which are bound to the membrane or the surface of a biomimetic object and therefore Y-junction 
formation is limited to the space close to the enzyme covered surface. 
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been employed previously to the actin network in the geometry of the lamellipodium 
ifnll . However, this model has a viscoelastic character and stresses build up only tran- 
siently due to temporary cross-link formation and not due to static crosslinks formed 
close to the polymerizing interface of the network. The rigorous derivation of a macro- 
scopic constitutive law of the network is therefore still lacking and is the subject of the 
present paper which takes up several ideas already outlined in 11811 . 

One crucial observation is that actin forms a more or less periodic network on the 
microscopic scale, which is stable on the time scale of symmetry-breaking (i.e. several 
minutes). The size of each elementary cell, e.g. the distance between two links in the 
actin network (~ several tens nm), is small compared to the total size of the structure 
(~1 um), which introduces a small parameter rj into the problem and allows for the 
upscaling of the microscopic mechanical model to a continuous medium by discrete 
homogenization methods 1 3 2(| 21 ] . 

In this paper we apply the discrete homogenization method to a simple 2D network 
structure which has been "grown" from a circle. While we keep the network structure 
and the microscopic mechanical model simple, we retain crucial biological features. 
In the Theory section we will define the network topology along with the microscopic 
mechanical properties and derive a set of continuous equilibrium equations. In the 
Results section we will solve these equations for the axisymmetric equilibrium state 
and identify the conditions under which linearly stable solutions exist. Finally, in the 
Conclusions section we will give some perspectives for this homogenization approach 
in modeling actin networks. 



2. Theory 

In this section we derive the continuous equilibrium equations governing the ho- 
mogenized network in an abstract Lagrangian coordinate frame. We start out with the 
discrete formulation of the network structure then upscale this description to a contin- 
uous medium. In a second step we introduce a free energy for the continuous medium 
and derive the equations of equilibrium using the first variation of the free energy. For 
completeness we outline in |Appendix A| and | Appendix B| the relation between the 
equilibrium description based on the stress vectors in the abstract Lagrangian frame as 
it is introduced in this Theory section and the classical Cauchy stress tensor. 

2.1. Description of the discrete network 

We consider a planar network around a solid circle of radius Rq. The network 
consists of stiff elastic filaments linked at their ends by nodes. Its topology is shown 
in Fig. |2] on the left. The bars are attached to the cylinder surface at N t sites evenly 
located at a distance p = r/Ro, i.e. at an angular distance r\ = 2n/N t . The net- 
work is made of N n layers of bars in the radial direction. As the network was formed 
by polymerization of actin monomers at the surface of the cylinder, each radial layer 
contains the same number of nodes N t . Each node can be identified by two integers 
(v x ,v 2 ) S {0, ...,N n - 1} x {0, ...,N t - l}as shown in Fig. |2] on the left. It is as- 
sumed that Nt and N n are very large and of the same order. Consequently, -q — 2ir /N t 
is assumed to be very small and the parameter a = ijN n is of order 1 with respect to 77. 
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Figure 2: [Color online] Sketch of a small part of the filament network with some examples for the node 
numbering before (left) and after (right) the transformation of variables (/j 1 ,^ 2 ) — > (A 1 , A 2 ) as intro- 
duced in Eqs. (6} and Q. The "curves" along which /i 1 =const. and /i 2 =const. are shown in blue and red, 
respectively. 



The elementary cell of such a network consists of one node and three bars b = 
1, 2, 3 and is shown in Fig. [3] Orientation and length of each bar is described by the 
vector B b . The tension in the bars is described by a linear law 



N" 



l b -l b 
l b 



(1) 



with l b = \\B b \\. k b is a spring constant and l b m is the length of the bar b at rest. Here 
we will only consider the situation k 1 = k 3 and l\ n = l 3 v At this stage the choice 
of a linear law is somewhat arbitrary and constitutes the simplest possible choice. It 
is consistent with basic principles in mechanics, such as material frame indifference. 
However, the homogenization procedure is valid for any microscopic law. 

So far we have only considered forces associated with the lengths of the filaments 
and have not included the fact, that we are only interested in networks which keep a 
certain topology, e.g. the elementary cell has a finite area || B 1 AB 3 ||> 0. To avoid this 
type of transition, i.e. a collapse of the network into zero area and a complete alignment 
of the vectors B b we add the following interaction moment between filaments 1 and 3 



with 



M 13 = M 13 B X A B 3 



ll 1 i i ) 2 (B 1 .B' 3 

M = k- 



:i\2 



(BKB 3 y 



(2) 



(3) 



where k is a constant and denotes the interaction strength. As will be seen later in 
the derivation of the equilibrium equations in section |2~3~1 expressions (O and ([3]) arise 
naturally from the first variation of the free energy which contains a potential which 
diverges for the complete alignment of B 1 and B 3 . 
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Figure 3: [Color online] Sketch of an elementary cell, showing the node numbering (u 1 , u 2 ) and the ele- 
mentary bar vectors B b (solid arrows). All other bars are shown as dashed lines. The "curves" along which 
i/^const. and i^ 2 =const. are shown in blue and red, respectively. 
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2.2. Homogenization 

The upscaling of the network to a continuous medium consists in determining the 
equivalent stresses from the bar tensions, the equations of equilibrium (or motion) sat- 
isfied by these stresses and an equivalent constitutive equation ensuing from the prop- 
erties of the bars. This can be carried out by using an asymptotic expansion (for an 
introduction see l22j|23l 12011 ). Here, as the network structure is simple, a more heuris- 
tic presentation can be used. The basic idea of the homogenization process is, that 
for most of the network motions, the positions of its nodes (y 1 , i/ 2 ) can be approxi- 
mated by a continuous deformation function ip (a 1 , /i 2 ) such that the position of the 
node (V 1 , v 2 ) is ip ((i lfl , /J 2?? ) with [i %r < = r\v % and (/i 1 * 7 , /i 2 '') <G {0, 77, . . . , a — 77} x 
{0, 77, . . . , 2tt — 77}. The discrete Lagrangian variables Ui 1,? , /i 2,) ) become then the 
set of Lagrangian (curvilinear) coordinates (fi 1 , /j 2 ) G lj of the equivalent continuous 
medium with u =]0, a[x]0, 2ir[. 

In this homogenized network the filament vectors B b can now be expressed using 
a simple Taylor expansion up to 0(rf) 



where d^i = -^j. Since they are associated with a quite simple numbering system 
for the nodes, the variables /i 1 and /i 2 arise naturally as Lagrangian variables of the 
equivalent continuous medium through the homogenization process. However, they are 
not the most convenient variables to study axisymmetric equilibrium configurations. 
Therefore we have introduced the variables A 1 and A 2 defined by 



and which correspond to the radial and angular position of the node in the network [see 
Fig.[2]to the right]. Setting 




(4) 



(5) 




(6) 




(7) 



one finds 




(8) 



with hi = d\itp and h 2 = o? where d\i = In the following we will employ 
the simplified notation di = d\i . 
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2.3. Equations of equilibrium 

The equations of equilibrium can be obtained directly by following the homog- 
enization procedure up to first order in r\ detailed in 12111 . However, here we have 
chosen a simpler approach using only the filament vectors B b obtained by homoge- 
nization and by defining a free elastic energy F of the network corresponding to the 
constitutive functions N b and the moment M 13 defined in ([T|i and (O. This free energy 
F can be expressed as the functional 



F[0\ = \ 



3 h.b 



(B 1 . B 3 ) 2 



.6=1 



2 (IH 3 ) 2 - (BKB 3 ) 2 



dX 1 dX 2 (9) 



where l b = \\B b \\. The first sum in (0 denotes the free energy associated with the 
compression and elongation of the filaments. The second term is associated with the 
interaction between the filaments 1 and 3 and serves as a "hard core" potential, which 
prevents the collapse of the networks, i.e. it prevents the complete alignment of the 
filaments 1 and 3. 

The first variation of F with respect to the deformation function reads 



SF 



1 



M 13 A 



dX 1 dX 2 
B 1 



3 R*> 

Y^N bi L.SB b + 



.6=1 
I .SB 1 



M 13 A 



- f dX 1 dX 2 ( S 1 .Shi + S 2 .Sh 2 



B 3 
(I 3 ) 2 



SB 



(10) 



where the stress vectors S 1 and S 2 are given by 



S 1 



N 1 
I 1 



-B 



N 3 
I 3 



■B 3 + M 13 A 



B l 



B 3 



N 1 iV 2 N 3 1 



B 3 



B L 



(I 1 ) 2 {I 3 ) 2 



(11) 
(12) 



Integrating Eq. (TlOb by parts and taking into account the periodicity in A 2 (i.e. the inte- 
grals along the boundaries of the network with A 2 = and A 2 = 2ir cancel each other 
out) one finds 



SF = - 



^ J (diS 1 + 9 2 5 2 ) .StpdX 1 dX 2 



i) 
1 

V Jo 



S 1 {a,X 2 ).SipdX 2 



V Jo 



S^O.X^.SipdX 2 



(13) 



The last two integrals over the boundary at A 1 =0 and A 1 = a in Eq. (TT3b evaluate to 
zero since tp is fixed at the internal boundary, i.e. Sip = for A 1 = 0, and there are 



7 



no normal forces acting on the external surface, i.e. an = S 1 /\\h 2 \\ = for A 1 = a. 
Here er denotes the Cauchy stress tensor (see |Appendix B i and ft is the unit outward 



surface normal vector in the deformed configuration. 

The mechanical equilibrium in the absence of volumic forces is therefore given by 

= SiS 1 + d 2 S 2 (14) 

with the boundary conditions as mentioned before 

VA 1 G ]0,a[, 0(X\O) = 0{X 1 ,2tt) 

VA 2 G ]0,2tt[, <^(0,A 2 ) = i? e r (A 2 ) (15) 
VA 2 G ]0,2tt[, ^(a,A 2 ) =0 

where e r = cos A 2 e\ + sin A 2 e*2 denotes the unit vector in the radial direction in the 
orthonormal coordinate basis {e*i, e^}. For a later use we introduce here also the tan- 
gential unit vector as eg = — sin A 2 e\ + cos A 2 e*2- 



3. Results 



3.1. Axisymmetric state 

In the following we will use the simplified notation tp' r = di(p r . Eq. ( fT4b with the 
boundary conditions ( fTBI ) was analyzed for axisymmetric solutions of the type (p = 
(p r (X 1 )e r - The stress vectors (fTTT i and (flZb can then be simplified to 



S 1 

S 2 = 7] 

The equilibrium equation ( TT4b becomes 



N 1 
I 1 
N 1 



V \ 2 — <p' r - M (Prtp'r \ e r 

N 2 



-<p r + M 13 <p r (p' 2 



= 8, 

with the boundary conditions 

<Pr(P) = i?0 



liV 1 



('0 



N 2 



(16) 
(17) 

(18) 



= 



N 



(19) 



X 1 =a 



Eq. (fT8l was solved using continuation techniques [24] and the solutions were charac- 
terized using the equilibrium network thickness H = (p r (a) — Rq and the associated 
free energy (9). 

At this stage we point out that the unknown function </? r (A 1 ) needs only to be 
piecewise continuously differentiable. It can have "corner points", i.e. discontinu- 
ities in the first order derivative. The presence of such discontinuous solutions can 



be inferred from the fact that d 2 f /dtp' 2 vanishes for a finite value cp' r at the "corner 
point" (A^, (fird)- For every continuously differentiable segment of an extremal solu- 
tion Eq. < f]~8T t holds supplemented by the following boundary conditions at the "free 
corner points" i25ll 

[df/d<ff r ] xl=x i_ = [df/d<p' r ]» =xi+0 (20) 
lf-<fr r df/d</ r ] xl=x i_ = If-dWdvQx^xi+o ( 21 ) 

The first condition < f20b is equivalent to the physical boundary condition, that the nor- 
mal stresses match at the "cornerpoint". The second condition (|2TT > minimizes the 
energy with respect to the position of the "corner point". 

First one might consider the special case k 2 = k = 0. As it is easy to see, iV 1 = 
solves Eqs. ( fT8b and ( fT9l . Under the condition that lj n > p/2 (recall that p denotes the 
distance at which the nodes are attached at the cylinder surface) one therefore finds the 
equilibrium solution 



R p 



sin 



A 1 ( p 



2 ' arCSm V 211 



(22) 



At the radial layer number A 1 = Aj = n — 2arcsin \j^~J one finds ip' r = and 
ifir/Ro = 21^/p, At this point the network collapses to zero volume and the filaments 
B 1 , B 2 and B 3 are oriented in parallel. The addition of further network layers will 
either lead to a back-folding of the gel into itself (in this case solution d22l holds), or 
the network layer will rest in this collapsed state, i.e. cp r /Ro = 2l} n /p for A 1 > A*. In 
this case A* constitutes a "cornerpoint", where Eqs. d20l ) and OTT) hold. 

The back-folding behavior is excluded in a physical situation since it requires the 
remodeling of the network structure and in the following we will not consider these 
special cases. The network collapse is unphysical as well since it implies a zero volume 
of the network layers. 

Now we will consider the more physical case k 2 > 0. Fig.|4]shows the dependence 
of the equilibrium network thickness H and the associated energy depending on the 
number of radial network layers a for the case k = k 2 /k 1 = 1 and = l 2 n = p, 
i.e. all filaments are identical and type 2 filaments are inserted at the internal interface 
at their equilibrium length. First, increasing the ratio of the elastic constant k 2 /k 1 
decreases the network thickness H for a constant number of radial network layers 
a. Second, the continuously differentiable solution branches with ip' r > cease to 
exist beyond a critical value a c . Instead solutions appear which have one or more 
discontinuities in tp' r and which contain solution segments with ip' r < (not shown). 

To exclude the back-folding behavior we have added an additional term in the free 
energy related to the interaction between B 1 and B 3 and which prevents the alignment 
of B 1 and B 3 . Fig. |4] shows exemplary the network thickness and free energy for a 
small value of n/k 1 . Solutions exist even for very large networks and the solution 
behavior for small networks is similar to the behavior with k = 0. 

Fig.[5]shows exemplary solutions for Eq. (TT~8T > for small, intermediate and large net- 
works. For small and intermediate networks, the filament lengths and stresses increase 
almost linearly with A 1 , whereas for large networks the behavior is strongly nonlinear. 



9 



0. 1 



0.8 1.2 

a 



1.6 



Figure 4: [Color online] (a) Equilibrium thickness of the network H and (b) the associated energy depending 
on the number of radial network layers a for k = and n/k l p = 0.025 as indicated in the legend. 
Remaining parameters are k2/k\ = 1 and = P- 



The network is always under radial compression, whereas it is under tangential com- 
pression close to the surface of the obstacle and under tangential extension at the outer 
surface of the network. 

3.2. Linear stability 

To study the (mechanical) linear stability of the axisymmetric state we have evalu- 
ated the sign of the second variation of F with respect to at the prestress axisymmet- 
ric state ip r e r 

5 2 F= [ Shi 4 ShjdX 1 dX 2 (23) 
Ju dhidhj 

with the boundary conditions 

VA 1 6 ]O,a[,60(X 1 ,Q) = 60(X 1 ,2Tr) 

VA 2 e ]0, 2tt[ , 6<p(0, A 2 ) = (24) 



Details of the stability analysis can be found in |Appendix C Briefly, we have used 
a Fourier mode decomposition of 60 and a direct integration in A 2 and finite elements 
for a numerical integration in the A 1 . 

Fig.|6]shows the linear stability of (p r (X 1 ) depending on the elastic ratio k 2 /k 1 and 
a. The linear stability of the axisymmetric state depends on the number of radial net- 
work layers superimposed over each other, whereby thin networks are stable and thick 
networks are unstable. Furthermore, networks with higher tangential stresses (large 
k 2 /k 1 ) destabilize at smaller a. The instability occurs independently of the transversal 
wavenumber m and is accompanied by the appearance of a region adjacent to the inner 
interface where d 2 f/dh 2 e < 0. The appearance of the instability is probably due to the 
fact, that as the network thickness increases the compressive load acting on type 1 and 
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Figure 5: [Color online] Solutions of Eq. ll8l for different network thicknesses a as indicated in the legend. 
Shown is (a) the radial node position, (b) the radial stress, and (c) the tangential stress distribution. Remaining 
parameters are k 2 /k 1 = 1, Kjk^p = 0.025, 1^ = = p. 



type 3 filaments in the innermost network layers increases. The elastic stress is released 
by rotating consecutive network layers against each other. Since the model does not 
contain any force opposed to a change in orientation of consecutive filaments of type 
1 and 3, which introduce an additional length scale related to a "persistence length" 
of consecutive filaments of the same type, a catastrophic growth of the instability at 
infinitely small scales occurs (data not shown). 



4. Conclusions 

We have introduced a two-dimensional homogenization model to describe the me- 
chanics of stiff elastic filament networks assembled on a cylindrical surface resembling 
Arp2/3 cross-linked actin networks as they occur in living cells. Our approach allows 
us to realistically model macroscopic deformations and stresses due to network growth 
by using only assumptions about the microscopic network properties. The main pur- 
pose of the model is to introduce a minimal set of mechanical ingredients as a basis for 
the more advanced problem of the time-dependent evolution of the network interfaces 
due to polymerization/depolymerization at the interface with a possible symmetry- 
breaking and subsequent motion. 

The approach is novel in the field of actin driven motion, inasmuch as it allows 
to calculate equilibrium stresses and the position of the interface between the network 



and the surrounding solution without any ad hoc assumptions |14J or restrictions to 
infinitesimal deformations Hill . Furthermore it will allow a rigorous calculation of the 
equilibrium state for networks with perturbed interfaces, resulting from spontaneous 
variations in the polymerization or depolymerization speed at the network interfaces as 



outlined in [18]. However, in this study we have limited ourselves to fixed interfaces 
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Figure 6: Linear stability of continuously differentiable axisymmetric solutions depending on the parameters 
k? /k 1 and a. Above (below) the solid line axisymmetric solutions are linearly unstable (stable). Remaining 
parameters are = = p and n/^p = 0.025. 

(in the Lagrangian frame) and studied the existence and linear stability of physical and 
biological relevant axisymmetric equilibrium states. 

We find that in a triangular network where only changes in the length of the edges 
contribute to the free energy the existence of sensible solutions is limited by the number 
of radial network layers a superimposed over each other. As one increases a beyond a 
critical number the network starts to fold back into the negative radial direction, since 
the extension of the microscopic filaments in the tangential direction is unfavorable 
as one moves away from the cylinder surface. The inclusion of an additional term in 
the free energy related to the interaction between two selected filaments prevented this 
back-folding. 

Furthermore, the linear stability of the axisymmetric solution depends on the mi- 
croscopic elastic properties. In general, thin networks are stable and thick networks 
are unstable, whereby a higher elastic constant for filaments oriented in the tangential 
direction favors the occurrence of the instability for thinner networks. The instability 
is independent of the transversal wavenumber of the harmonic perturbations. A closer 
inspection of the instability reveals that it occurs first at the cylinder surface and is 
reminiscent of a buckling instability of type 1 and 3 filament, i.e. a counterrotation 
of consecutive radial network layers. However, there is no mechanism in this simple 
model to stabilize this instability. A more advanced approach would have to include a 
more detailed microscopic mechanical model with nonlocal contributions (i.e. in the 
spirit of second gradient model), which is beyond the scope of this paper and will be 
the subject of further investigation. 
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Figure A. 7: A facet S with the unit normal it and the segment length s is placed in the network. The 
elementary cell with the basis vectors B b and the height s is shown in gray. 



Appendix A. Derivation of the Cauchy stress tensor in the network 

The first to define a stress tensor in a discrete medium in the Lagrangian frame were 



Cauchy and Poisson whose derivation can be found in 123112211 . We present here the 
basic ideas in two dimensions to facilitate the comprehension of the Theory section. 

For the derivation of the stress one considers a periodic structure of nodes and 
their connections. Forces act only along the directions of the node connections. The 
periodic structure is divided into elementary cells which are numbered by (V 1 , v 2 ). The 
connections are identified by the triplet (b, v 1 , v 2 ), where b identifies the connection 
within the elementary cell. 

The stress tensor is defined over a force R acting on a arbitrary facet S with unit 
normal n placed in the periodic medium as shown in Fig. |A.7| The force R is obtained 
by summing up the forces acting along the node connections which are cut by the facet. 
Thereby one assumes that the length S of the facet is sufficiently small that the forces 
exerted by identical connections are constant over the facet. Therefore the force R is 
given by 

B b 

R = V mhga(B b .n)N h . (A.l) 
V \\ Bb W 

Here N b and m denote the tension and number of node connections of type b. The 
term sgn(B b .n) accounts for the orientation of the connections. There are different 



ways to determine m h , we give here only the result derived in detail in 112211 . The basic 
idea is, that the facet S is cut by the vectors of one type b into segments s and therefore 
m b = S/s. As can be seen in Fig. IA. 71 for the example of the type 1 filament, the 
height of the elementary cell is s = slB 1 .ril/WB 1 ]] and consequently its area is given 
by A c — sH^B 1 1| . Therefore the number of links of type b cut by the facet S can be 
expressed as 



m b = ~ i -i (A _ 2 ) 
Ac 



S\B 



b 
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Using expression dA.2b with Eq. ( IA. 1 b one obtains 

R= -J- Y N b -^-(B b .n)S . (A3) 

\\ Bb \\ 

After some rearrangements and using the identity 

R = anS (A.4) 

one finds 

IW*^ (A.5) 
A c ^ \\B"\\ 

Appendix B. Relation between the stress vectors S x and the Cauchy stress tensor 

In the present section we demonstrate the relation between the classical equilib- 
rium equation known from continuum mechanics using the Cauchy stress tensor er in 
the deformed configuration and the equilibrium formulation we have used throughout 
this paper defined on the Lagrangian reference configuration. Note, that we will use 
the Einstein notation, i.e. summation over repeated indices, and that we will use the 
abbreviation dj = d\* which denotes a partial derivation with respect to the coordinate 
X 1 in the Lagrangian reference configuration. 

We start out by introducing some equalities. Let <p be the deformation mapping of 
the network with respect to its Lagrangian configuration uj and let x = tp (A 1 , A 2 ) be 
the position in the deformed configuration ui^ = (uj) of the network point (A 1 , A 2 ) 
of the Lagrangian reference configuration. Let dx be the difference of positions of two 
neighbouring points (A 1 , A 2 ) and (A 1 + dX 1 , A 2 + <iA 2 ), then 

dx = h l dX l (B.l) 

where hi = di<p. Consequently one finds 

dX l =h\dx (B.2) 

with the property 

hi.h j =5i (B.3) 

Let now v be a virtual velocity field which, due to the one to one mapping Cp from 
u to ui v , can be considered either as a function of A 1 , A 2 or as a function of x. v is 
assumed to be differentable and let Vu be its differential with respect to x and <9,i7be 
its partial derivative with respect to A 1 , then 

dv = Vvdx (B.4) 
dv = d t vdX l (B.5) 
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Now, the use of equations ( IB. lb and (lB.2t yields 

Vv = d t v®h l (B.6) 
diV = Wvhi (B.7) 

In the absence of body forces the virtual power formulation of the equilibrium 
equation ( fPfl i reads 



S i .d i vdX 1 d\ 2 = (B.8) 

where the virtual velocity field v vanishes on the boundaries A 1 = and A 1 = a of w 
and is 2ir periodic with respect to A 2 . 

The change of variables (A 1 , A 2 ) ■(-> x = (^(A 1 , A 2 ) in the intergral of the previous 
equation yields 

S\{Vvhi) -ds = (B.9) 

— <p 

where ds — h dX 1 d\ 2 denotes the area element of the body in the deformed configu- 
ration with h = \\hi A h%\\. The previous equation also reads 

j i (& ® : Vvds = (B.10) 

where A : B denotes the usual twice contracted product between the second order 
tensors A and B. 

Comparing with the classical virtual power formulation of the equilibrium in the 
Eulerian variable x which reads 

a:Vvds = (B.ll) 



we can see that the expression of the Cauchy stres tensor er in terms of the S^'s reads 

a = \s l ®h l (B.12) 
h 

Conversely, the use of definition (IB. 31 of the contravariant base (ji 1 , h?^j in the 

equation (IB. 12b yields the expressions of the vectors S 1 and S 2 in terms of er which 
reads 

S l = hah 1 (B.13) 
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Appendix C. Linear stability analysis 

Since the boundary conditions are periodic in A 2 one can decompose 50 = 5(p r e r - 
Sipeeg into a Fourier series 



5<p r = [^™(A 1 )cos(mA 2 ) + ^™(A 1 )sin(mA 2 )] 

m=0 

oo 

5<p e = ^ [^(A 1 )cos(mA 2 ) + <5^(A 1 )sin(mA 2 )] 



(CI) 



where m denotes the wavenumber and <5</?™, Sip™, 5tp™ c , Sip^\ denote perturbation am- 
plitudes which depend only on A 1 . Using ansatz dC.U one can perform the integration 
of (f23) in A 2 directly whereby different wavenumbers m decouple. One finds then 

OO „Q, 

S 2 F = 7rV / (S 2 f? + 5 2 f 2 n ) dX 1 (C.2) 

where 5 2 f{ n 2 takes the form 

6 2 f™ 2 = A(a') 2 + B{b') 2 + 2C(a'a + ma , b)- (C.3) 
2D(b'b + mab') + E(ma + b) 2 + F(a + mb) 2 

with (a, b) = (Stp^Sip™) and (a, 6) = (-5<p™,5<p%,) for <5 2 /f and d 2 f^, respec- 
tively. The factors A, . . . , F are evaluated at the axisymmetric solution = ip r e r and 
are given by 



A dh\ r - p +Zk (i^r +K tfrf 

B 92 f -2 Nl I ^ 2,: ^-^ 

dh\ B 2{l i )3 ^ 

d 2 f _ , tfrrt ^/16 + < 4 

dh lr dh 2 g 2(P)3 iffy? 

n d 2 f lVrV ' r (y 2 /4-^ 2 ) 2 

dhxedhir 2(Z!) 3 + ipltp? 

d*f _ iV 1 , AT 2 , , x y; 2 , ^/4^; 2 

2Z 1 + / 2 + 2(^) 3 + 2^; 2 



9^ 2 e 2Z 1 Z 2 8(^) 3 (Z 2 ) 3 



K - 
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The integral iC.2i was calculated numerically by using finite elements, whereby the 
functions 8(p™, 6cp™ c , Stp™ s , and Sip^, were approximated by Pi-elements and the func- 
tions A, . . . ,F were approximated by Po elements. This results in the quadratic form 

oo 

* = 2^ XVPrcVes)^ Wrofffa) +{-^3^00)0, {-frs^ec) J ( C - 4 ) 
m=0 

where " indicates a vector containing function values evaluated at the discretization 
points A|, (i = 1, 2, . . . , n) and Q m is a symmetric 2n x 2n matrix. 

For the axisymmetric state to be linearly stable S 2 F > 0, i.e. Q m needs to be 
a positive definite matrix. We have tested this condition for wavenumbers m < 4. 
But it turned out that an instability appears independently of m when the factor B = 
d 2 f/dhf e changes sign at A 1 = 0. The corresponding phase diagram is shown in 
Fig.E] 
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